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1. Introduction 

The Standard Model (SM), by meeting all confrontations with experiments, stands 
as a remarkably simple parameterisation of known physics. Yet it has many un- 
satisfactory aspects such as a large number of parameters, a triplication of chiral 
families, and three different gauge structures. Consequently, it is strongly believed 
among theorists that there must exist a simpler underlying structure of which the 
SM is the low energy piece. Many think that such a structure will make its appear- 
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ance at much higher energies, somewhere in the unexplored region between a few 
TeV and the Planck scale. Experimentalists can only proceed one or two orders of 
magnitude at a time in their exploration of those scales. Theorists, on the other 
hand, have mainly a twofold approach in divining this structure. In the first, the 
quantum numbers of the SM are grouped into mathematically pleasing structures 
(constrained however by the predictions affecting the low energy part of the theory), 
resulting in an exercise in quantum pattern recognition. This has led to Grand Uni- 
fied Theories (GUTs) pQ. The second approach is to use the renormalisation group 
to extrapolate the SM parameters to the unexplored scales [2]. The purpose is to 
find if those parameters satisfy interesting relations at higher energies. When used 
in conjunction with the former approach, this can give powerful hints of the physics 
expected at higher energies. 

We shall also recall that in order to obtain finite results in quantum field the- 
ory, in a higher order than the tree level, one has to perform the renormalisation 
program. The independence of the renormalisation procedure from the renormali- 
sation point leads to the dependence of the Lagrangian parameters on the point of 
renormalisation. This dependence is governed by the Renormalisation Group Equa- 
tions (RGEs) for the coupling constants and other parameters of the Lagrangian. 
As such, it all depends on having accurate data to input as initial conditions on the 
RGEs, as well as a strong theoretical basis for the evolution equations themselves. 
These theories, at asymptotic energies, may reveal some new symmetries or other 
interesting properties that give deeper insight into the physical content. Also, from 
the requirement of the stability of the theory one can, for example, give bounds on 
the physical parameters, like the Higgs mass in the SM or its extensions. 

In the SM, the runnings of the gauge, Yukawa and quartic scalar couplings is 
logarithmic with the energy scale. Although the gauge couplings do not all meet 
at a point, they tend to unify near 10 15 GeV. Extensions to the SM such as extra- 
dimensional scenarios accessible to SM fields have the virtue, thanks to the couplings 
now having a power law running, of bringing the unification scale down to an ex- 
plorable range [3] . Note that many other extensions to the SM exist which alter the 
runnings in different ways, such as supersymmetry, where a range of new particles 
ensure the gauge couplings do meet at a point, but runnings remain logarithmic. 
However, combinations of the different approaches have led, in the last two decades, 
to some extremely interesting phenomenological models, such as TeV scale extra- 
dimensional scenarios giving rise to new supersymmetry breaking mechanisms [4], 
as well as addressing fermion mass hierarchy [5] and providing cosmologically viable 
dark matter candidates [6]. 

With the Large Hadron Collider (LHC) now up and running, exploration of the 
realm of new physics that may operate at the TeV scale has begun. Among these 
models those with extra spatial dimensions might be revealed in such higher energy 
collider experiments, where the Universal Extra Dimension (UED) model makes for 
an interesting TeV scale physics scenario; as it features a tower of Kaluza-Klein (KK) 
states for each of the SM fields, all of which have full access to the extended space- 
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time manifold [7 . This particular scenario has recently been extensively studied 
in the literature, such as investigations of electroweak symmetry breaking, proton 
stability, gauge hierarchy and fermion mass hierarchy problems, B physics, dark 
matter etc. [&Ml3|. This model has been a fruitful playground for addressing a variety 
of puzzles in the SM. 

We therefore collect in a comprehensive manner and in one place the necessary 
tools for making renormalisation group analyses of the SM and its UED extensions. 
We view this compendium as a template for general applications of the running of 
parameters from the mz scale to the Planck scale. The observable parameters of the 
SM are: 6 quark masses, 3 lepton masses, 4 parameters of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [14] and 3 gauge couplings. The RGEs for the CKM matrix 
being obtained from the RGEs for the Yukawa couplings. This can also be extended 
to include neutrino masses and mixings possible in the leptonic sector, as we shall 
discuss further in section |9] 

In this review we first introduce the various models and their varieties we shall 
consider (section[2]and section[3]for the supersymmetric extensions to this), next de- 
riving the RGEs for the gauge couplings constants (section |4} and Yukawa couplings 
for the SM and UED scenarios (section [5]) and 5-dimensional Minimal Supersym- 
metric SMs (5D MSSM) (section EJ. This shall be followed by a review of the CKM 
parameters evolution (sections [7] and [8|). Extensions to massive neutrino scenar- 
ios and their mixings evolution will follow in section [9] With a summary and our 
conclusions in section [TOl 

2. The UED Standard Model 

The UED model places particles of the SM in the bulk of one or more compactified 
extra dimensions. In our case we have a single flat extra dimension of size i?, com- 
pactified on an S1/Z2 orbifold. As such we will have an infinite tower of KK modes 
with the zero modes corresponding to the SM states. These KK modes are in the 
TeV scale and modify the running of the RGEs at relatively low energy scales. The 
UED model, like any higher dimensional theory, is an effective field theory which 
is valid up to some scale A, at which a new physics theory emerges. Between the 
scale where the first KK states are excited and the cutoff scale A, there are 
finite quantum corrections to the Yukawa and gauge couplings from the AR num- 
ber of KK states. Up to the scale the first step KK excitation occurs, the RG 
evolution is logarithmic, controlled by the SM beta functions. With the increasing 
of the energy, that is, when the KK threshold is crossed for each successive mode, 
new excitations come into play and govern new sets of beta functions. The values 
of physical parameters such as Yukawa couplings and gauge couplings do not run in 
the old SM fashion, instead they receive finite quantum corrections whose magni- 
tudes depend explicitly on the value of this cutoff parameter. As a result, once the 
KK states are excited, these couplings exhibit power law dependencies on A. This 
can be illustrated if AR ^> 1, to a very good accuracy, the generic SM beta function 
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is shown to have the power law evolution behaviour |10j : 



(1) 



where (3 is a generic contribution from a single KK level, and where its coefficient 
is not a constant but instead S(fi) = fiR, with n Max = A, reflecting the power law 
running behaviour. As a result of faster running, the gauge couplings tend to lower 
the unification scale down to a relatively low order, which might be accessible to 
collider experiments, such as the LHC or the proposed International Linear Collider 
(LLC). Therefore, constraints from precision electroweak tests and current (or fu- 
ture) collider data would yield bounds on the compactification radius R. The RGE 
are an important tool for the search of the properties of the quark masses and the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix at different energy scales. It is there- 
fore of great interest to have an implementation of the UED model in studying these 
RGE. 

The first version of this model we shall consider, the bulk UED model, sometimes 
known as the minimal UED model, has one extra dimension compactified on a circle 
of radius R with a Zi orbifolding which identifies the fifth coordinate y — >■ — y. The 5- 
dimensional KK expansions of the weak doublet and singlet as well as the Higgs and 
gauge fields are shown (the corresponding coupling constants among the KK modes 
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The zero modes in the above equations are identified with the 4-dimensional SM 
fields, whilst the complex scalar field H and the gauge field are Z2 even fields, 
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and there is a left-handed and a right-handed KK mode for each SM chiral fermion. 
Note that in models with UED momentum conservation in the extra dimensions, we 
are led to the conservation of KK number at each vertex in the interactions of the 
4-dimensional effective theory (or strictly speaking, the KK parity (—1)" is what 
remains conserved, where n is the KK number). In the bulk we have the fermion 
and gauge field interactions as follows: 

7T-R 

C-Leptons = J dy{iL(x, y)T M V M L(x, y) + ie(x, y)r M V M e{x, y)} , 
o 

ttR 

^Quarks = J dy{iQ(x, y)T Mr D M Q(x, y) + iu(x, y)T M V M u{x, y) + id(x, y)T M V M d(x, y)} , 

(3) 



i) 



where T M = (7^, 17 s ), and M = 0, 1, 2, 3, 5. Explicitly, the kinetic terms are given 
by: 

T> M Q(x, y) = (d M + iglG M + ig\w M + i^glBu) Q(x, y) , 

V M u(x, y) = (d M + ig\G M + *^.9i^a/^ u(x, y) , 

T> M d(x, y) = (d M + iglG M + i-^-g\BM^ d(x, y) , 

V M L{x, y) = (d M + ig\W M + i^glBjA L(x, y) , 

T> M e(x, y) = (d M - ig\B M ) e(x, y) . (4) 

The gauge couplings gf, g\ and g\ refer to those of the SU(3), SU{2) and U(l) 
gauge groups respectively, and are related to the 4-dimensional SM coupling con- 

g 5 

stants by gi — 1 and the five dimensional gauge fields have the generic 
V ttR 

form Am = (A^^Ac,). After integrating out the compactified dimension, the 4- 
dimensional effective Lagrangian has interactions involving the zero mode and the 
KK modes. However, these KK modes cannot affect electroweak processes at tree 
level, and only contribute to higher order electroweak processes. For the one-loop 
diagrams of the Yukawa couplings we choose the Landau gauge in what follows, as 
many one-loop diagrams are finite in the Landau gauge and have no contribution 
to the renormalisation of the Yukawa couplings. We therefore consider the RGE 
for the quark- Higgs Yukawa couplings from which we obtain the evolution of the 
quark masses and the CKM matrix. The one-loop Feynman diagram contributions 
to the Yukawa couplings in the SM and UED model have been explicitly illustrated 
in [T0lll5fH7] . In the UED model, where for each energy level rii, we effectively have 
a heavier duplicate copy of the entire SM particle content. That is, at each KK 
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excited level, the KK tower corresponding to the fields in Eq. @ exactly mirror the 
SM field ground states. However, new contributions from the A$, 



interactions (that of the fifth component of the vector fields, i.e. A$ = {G5, W5, -B5}) 
also contribute. In contrast, the fifth component of the gauge bosons A§(x,y) is a 
real scalar and does not have any zero mode, transforming in the adjoint represen- 
tation of the gauge group. 

A simple alternative to this model is that of the brane localised UED model, 
where we have the same fields but where the fermion matter fields cannot propagate 
in the bulk and they are restricted to the brane. In the brane localised UED model, 
extra dimensional models make an interesting TeV scale physics scenario and as 
discussed earlier might be revealed in higher energy collider experiments, as they 
feature states that have full access to the extended space-time manifold. For the 
case of brane localised matter fields, only the boson fields (the gauge fields and 
the scalar fields) can propagate in the bulk space. However, if the compactification 
radius R is sufficiently large, due to the power law running of the gauge couplings, 
it will enable us to bring the unification scale down to an exportable range at the 
LHC scale. 

3. The 5D MSSM 

Another useful model we shall consider is the 5D MSSM defined in [31 HM26] . The 
5D MSSM is a five dimensional N — 1 supersymmetric model compactified on the 
S1/Z2 orbifold which breaks the 5D Lorentz invariance to the usual 4D one. This 
breaking gives a momentum conservation along the fifth dimension which conserves 
the KK number at tree level and KK parity at loop level. One of the main implica- 
tions of KK-parity invariance is that the lightest KK mode is stable and can be a cold 
dark matter candidate. In this compactification we can recover the MSSM at zero 
mode since we obtain chiral fermions. Note though, that since we are working with 
a supersymmetric model, the beta function can be derived in the superfield formal- 
ism, where we shall discuss Af = 1 supersymmetry in a five-dimensional Minkowski 
space and its description in terms of 4D superfields. The space-time coordinates 
being labelled by (x^,y). 

The gauge sector is then described by a 5D J\f = 1 vector supermultiplet which 
consists (on-shell) of a 5D vector field A M , a real scalar S and two gauginos, A and 
A'. The action for which can be given by: 




71=1 




Sg 




i\'r M D M X' + (A + A') [5, A + A' 



(6) 
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with D M = d M +iA M and T M = (7^, ij 5 ). F MN = -±[D M ,D N ] and k normalises 
the trace over the generators of the gauge groups. 

From the decomposition of the 5D supercharge (which is a Dirac spinor) into two 
Majorana-type supercharges, which constitute a Af = 2 superalgebra in 4D, one can 
rearrange these fields in terms of a Af — 2, 4D vector supermultiplct, O = (V ,%): 

• V : Af = 1 vector supermultiplet containing A M and A, 

• x '■ N = 1 chiral supermultiplet containing A' and S' = S + iA 5 . 

Both V and \ (and their component fields) are in the adjoint representation of the 
gauge group Q. Using the supermultiplets one can write the original 5D Af = 1 
supersymmetric action Eq.© in terms of Af = 1 4D superfields and the covariant 
derivative in the y direction [21j : 




with W a = ~\D er 2gV D a e 2gV . D a is the covariant derivative in the 4D TV = 1 
superspace (see Refs [27ll28] .l and V y = d y +\- To find the Feynman rules to a given 
order in the gauge coupling g, one can expand and quantise the action [18]. The 
beta functions for the couplings of the operators in the superpotential are governed 
by the wave function renormalisation constants Zij = 1 + SZij due to the non- 
rcnormalisation theorem [29j . The Feynman diagrams related to the wave-function 
renormalisation are given in FigfT] 

The Higgs superfields and gauge superfields will always propagate into the fifth 
dimension. Different possibilities for the matter superfields will be discussed, where 
superfields containing SM fcrmions can propagate in the bulk or are restricted to 



a) n (0) b) n c) (0) 




n (0, n) n (0) 



d) (0) e) m 




Fig. 1. The one-loop diagrams related to the wave-function renormalisation of the matter super- 
fields, in which diagrams a)-e) refer to the case where all the matter fields are in the bulk, and the 
excited KK states are labelled by the number without the bracket; whereas diagrams a), c) and 
d) are related to the brane localised matter fields case, in which the KK states are labelled by the 
number inside the bracket [181 1231 . 
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the brane. For the case where all fields can propagate in the bulk, the action for the 
matter fields would be [18] : 

d 8 zdy + <F=$ 4 C + ^8^5(6) - ^8^5(0) 

+g(2^V^-2^V^ + ^X^S(0) + ^x^(0))} ■ (8) 

Again, this action can be expanded and quantised. The x-field should be odd under 
Z2 symmetry because it appears together with a derivative d y , whereas V is even. 
For the two matter superfields, we choose $ to be even and the conjugate $ c to be 
odd such that $ c vanishes on the brane. Only the even fields have zero modes. The 
Fourier decomposition of the fields being: 
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We can write the action for the second case where all superfields containing SM 
fermions are restricted to the brane. In which case the part of the action involving 
only gauge and Higgs fields is not modified, whereas the action for the superfields 
containing the SM fermions becomes: 



S n 



= J d 8 zdy6(y) + 20iV$i} 



(10) 



Due to the 5D M = 1 supersymmetry, Yukawa couplings are forbidden in the 
bulk. However, they can be introduced on the branes, which are 4D subspaces with 
reduced supersymmetry One can also add the effective neutrino mass operator 
(also called lepton number violating Weinberg operator) , with dimensional coupling 
kij in which we are interested to show its evolution and therefore the Majorana 
mass term for neutrinos. We will write the following interaction terms, called brane 
interactions, containing Yukawa-type couplings: 



Sbrane = / ^ zdy8(y) 



^A,,;."!\<!yh, 



Lj Hit L i Hn 

AM 



6(0) + h.c. 



(11) 



where L and H u are the lepton and up-type Higgs doublet chiral superfields respec- 
tively. This operator is used to study neutrino masses and mixings, where RGEs 
for this effective operator have been derived in the context of the four-dimensional 
SM [30] and MSSM [31] and shall be discussed further in section [9] An extension 
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to compactified extra-dimensions was considered in Ref |18) . , and we shall use a 
similar formalism in the next section. 

From FigQ]we obtain the wave function renormalisation of the matter superfields 
and then we calculate the beta function. If we exclude the effects of extra dimensions, 
one can obtain the beta function of the usual MSSM. The Figs. QJO-e) refer to 
the case where matter superfields are able to propagate in the bulk (excited KK 
states are labelled by the number without brackets) and diagrams in FigsQ]a), c )> 
d) represent the case of brane localised matter superfields (excited KK states are 
labelled by the number inside brackets). 



4. Gauge couplings 

The evolution of the gauge couplings in four dimension at one loop are given by: 

16tt 2 ^ = b l9 * , (12) 

where bf M = {%, -f , -7) and bf ISSM = (f , 1, -3) [32], using a SU(5) normalisa- 
tion. If we consider our 5D theory as effective up to a scale A, we have contributions 
from the KK modes which give a power law evolution to the gauge couplings since 
extra-dimensions cause their running to vary much more rapidly. Eq.(JT])can be writ- 
ten in terms of the scale parameter t: 

16vr 2 % = [b % + (S(t) - l)k]g t 3 , (13) 
at 

where bi take the following form in the case of the model UED SM [151133] : 
,7 r 7 x ( 1 41 21\ 8 

(b 1 ,b 2 ,b 3 )=( ro ,- 1 r,- Y )+- 3 V, d4) 

with r] being the number of generations of matter fields in the bulk. Therefore, for 
all our matter fields propagating in the bulk (that is, 77 = 3), we have for the UED 
SM bulk case: 

'81 7 5 N 
10' 6' _ 2, 



(15) 



Similarly, for all our matter fields localised to the 3-brane (that is, i] — 0), we have 
for the UED SM brane case: 

/ 1 41 21\ 

^(lo'-T'-Tj ' (16) 

Next we consider the beta functions of the gauge couplings in the 5D MSSM. In 
fact, after compactification of the 5D MSSM, we have two 4D Af = 1 chiral super- 
multiplets, $ and $ c , where the zero modes of $ give us the normal matter fields 
and Higgs fields as well as their super partners, while $ c is a new supermultiplet. 
In the component field formalism, at each KK level, aside from the quantum cor- 
rections that mirror those of the 4D MSSM, the only new one-loop contributions to 
the Feynman diagrams are from the wave function renormalisation of (which 
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contribute via the coupling of to the complex scalar field and its super-partner 
in the superfield \, and the coupling of with the new fermion field and its super- 
partner in the superfield $ c associated with the two doublets of the Higgs fields and 
the matter fields respectively in the bulk). This then gives rise to the master beta 
functions of the gauge couplings in the 5D MSSM as follows [23] : 




where rj represents the number of generations of fermions which propagate in the 
bulk. So in the two cases we shall consider, that of all fields propagating in the bulk 
(rj = 3) we have [3~4] : 



Similarly, for all our matter fields localised to the 3-brane (that is, rj = 0), we have: 



In Figs [2] and [3] we have plotted the running of the gauge couplings for the 
UED SM case and the 5D MSSM respectively for the brane localised and bulk field 
cases, and for several choices of compactification scales for the extra- dimension 
(R). From these plots, and the discussion given in Ref [33]. , we find that for 
the three gauge coupling constants to approach a small region at some value of 
t requires an extremely large value of 1/R, which is of no phenomenological interest 
at present. For the case of our fields being brane localised in the UED model, 
we see a similar behaviour: the extra-dimensions naturally lead to gauge coupling 
unification at an intermediate mass scale for the compactification radii considered 
here. Furthermore, as illustrated in Figj2] the extra space-time dimensions naturally 
lead to the appearance of GUTs at scales substantially below the usual GUT scale. 

We assume the fundamental scale is not far from the range of LHC scale, and set 
the compactification radii to be R^ 1 = 2 TeV, 8 TeV, and 15 TeV respectively. In the 
limit when the energy scale is much smaller than R~~ 1 , since the energy of the system 
is less than the excitations of the first KK modes, the theory reduces to the usual 
4-dimensional SM, and the existence of the KK excitations are ignored. When /i > 
excitations of many KK modes become possible, and the contributions of these 
KK states must be included in all physical calculations. This is characterised by the 
second term in Eq.([T]) in the general beta function. Once the energy passes R^ 1 the 
excited KK modes tend to increase rapidly the running of the gauge couplings, and 
ultimately change the scale dependence of the gauge couplings from logarithmic to 
those of a power law as a function of /i. Quantitatively, due to the fast running of 
the gauge couplings, we find they nearly meet at around t = 6.4, 7.8, 8.4 for radii 
R^ 1 = 2, 8, 15 TeV respectively. The extra dimensions naturally lead to gauge 
coupling unification at an intermediate mass scale. 





(17) 




(18) 




(19) 
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Similarly, in Figj3]we have plotted the evolutions of the brane localised and bulk 
field cases for several choices of compactification scales for the 5D MSSM. From 
these plots, and as noted above, we find that for the three gauge coupling constants 
to approach a small region at some value t requires an extremely large value of 
l/i?, whereas, for the case of our fields being brane localised, the extra dimensions 
naturally lead to gauge coupling unification at a similarly valued intermediate mass 
scale for the compactification radii considered here. 

5. Yukawa evolutions in the SM and UED SM 

In the quark sector of the SM, we have ten experimentally measurable parame- 
ters, i.e. six quark masses, three mixing angles, and one phase (these angles and 
phase being encoded in the CKM matrix which we shall discuss in section [7]). A 
completely satisfactory theory of fermion masses and the related problem of mixing 




3456789 10 3456789 10 



Fig. 2. Gauge couplings (gi (red), 32 (blue), 33 (green) with: in the left panel, all matter fields 
in the bulk (UED bulk); and the right panel for all matter fields on the brane (UED brane); for 
three different values of the compactification scales (2 TeV (solid line), 8 TeV (dot-dashed line), 
15 TeV (dashed line)) as a function of the scale parameter t in the UED SM. 




Fig. 3. Gauge couplings (gi (red), 92 (blue), 33 (green) with: in the left panel, all matter fields 
in the bulk; and the right panel for all matter fields on the brane; for three different values of the 
compactification scales (2 TeV (solid line), 8 TeV (dot-dashed line), 15 TeV (dashed line)) as a 
function of the scale parameter t in the 5D MSSM. 
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angles is certainly lacking at present, however, there has been considerable effort to 
understand the hierarchies of these mixing angles and fermion masses in terms of 
the RGEs [171132 , 35 - 39j . First though we must recall that in order to explore the 
physics at a high energy scale we use RGEs as a probe to study the momentum 
dependence of the Yukawa couplings, the gauge couplings, and from these the CKM 
matrix elements themselves. As such we can consider one of the primary goals of the 
LHC as being to uncover any new dynamics within the TeV range, where instead 
of assuming the RGE goes from the Mz scale up to the GUT scale (10 15 GeV) by 
using the SUc(3) x SUl(2) x UyiX) symmetry, we know that models with extra 
dimensions may bring down the unification to a much lower energy scale. However, 
when using the RGEs as a probe, the initial values we shall adopt are very impor- 
tant, where we shall scale for the gauge couplings and the fermion masses at the 
Mz scale are shown in Table [TJ 



Table 1. Initial values for the gauge couplings, fermion 
masses and CKM parameters at Mz scale. Data is taken 
from Ref ISgPOlHI] . . 



Parameter 


Value 


Parameter 


Value 


at 


0.01696 


m e 


0.48657 MeV 


a 2 


0.03377 


m M 


102.718 MeV 


OL3 


0.1184 


m T 


1746.24 MeV 


m u 


1.27 MeV 


\v ub \ 


0.00347 


m c 


0.619 GeV 


\v cb \ 


0.0410 


mt 


171.7 GeV 


\Vus\ 


0.2253 


m d 


2.90 MeV 


J 


2.91 x 10 -5 


m s 


55 MeV 






mb 


2.89 GeV 







Furthermore, we shall also attempt, in section[9j to develop the RGEs of the lep- 
ton sector (including possible mixing angles and phases), which will require knowl- 
edge of the evolution of a parameter k, where the lowest order operator which gen- 
erates Majorana neutrino masses after electroweak symmetry breaking (EWSB), 
is the lepton-number violating Weinberg operator |42j . This lowest order operator 
(appearing with dimension d = 5 in four space-time dimensions) can be written as: 

-^(L c Je^^)(Lie^^) + h.c, (20) 

where L and 4> are the lepton and the Higgs doublet fields. M is the typical heavy 
energy scale for the range of validity of the low-energy effective theory. An operator 
of this type can be generated, for instance, by the usual see-saw mechanism. In 
which case the heavy scale M can be identified with the mass of the heavy right- 
handed neutrino. After EWSB the Higgs acquires a vacuum expectation value (vev) 
and the operator in Eq. (|20[) gives a Majorana mass term for the neutrinos. In the 
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context of the MSSM it can be written in the form: 

- ^(L^H^(L^H^) , (21) 

where L and H u are the lepton and up-type Higgs doublet chiral superfields respec- 
tively. This operator is crucial for the study of neutrino masses and mixings, where 
renormalisation group equations for this effective operator have been derived in the 
context of the four-dimensional SM and MSSM [301131] . 

In the present case we consider the effective neutrino mass operator with dimen- 
sional coupling kij] after spontaneous symmetry breaking, the Majorana neutrino 
masses can be written as m v = kv 2 sin 2 /3 (v being the vev of the Higgs field and 
tan/3, the ratio of the vevs of our two Higgs doublets) and k — kij/(2MwR) for 
bulk propagating, and k = kj (2M) for brane localised matter superfield scenarios 
respectively. 

As such, we have set Mz as the renormalisation point, and use t = ln( ■£-) and 
S(t) = e t AIzR- The general form of evolution equations for Yukawa couplings and 
neutrino k coupling at the one loop can be written in the following form Refs [43H45] 



16tt 2 ^ - Y d |r t Ci -G d + ^(YjY d - F„ty„)C 2 | , 
16tt 2 ^ = Y u jllCi - G u + ^{Y^Y U - yjr d )C 2 J , (22) 
16tt 2 ^ = Y e (T t d - G e + |(F e ty e )C 2 | , 
ak + {[Y e T Y e *]k + k[YjY e ]) C 3 . 



16tt 2 ^ 



, dk 



where T* = Tr 



3YjY d + 3YjY u + Y^Y e 



5.1. Standard Model 

The SM is a limiting case for the UED, where the KK states decouple. The coeffi- 
cients in the evolution equation are defined by: 



(23) 



GdSM 


= (\gi + \gl + 8<?I 


GuSM 


A7 2 9 2 
= U^ + 4 52+8 - 9 


GeSM 


2 9 2 \ 

= [jgi + i 9 2 J , 


asM 


= 2 T t - 3g 2 + A , 


GlSM 


= 1 , 


G2SM 


= 1 , 


G^SM 


3 

~~ ~2 ' 



11 



Note that there is a difference in the coefficients of the gauge coupling g\ between 
Ref [43] . and Ref [44]. due to the SU(5) normalisation factor 3/5. 

5.2. UED SM Bulk 

The UED contribution is obtained when KK states enter, where due to the orbifold- 
ing the zero mode for fermions are chiral, which are replaced by Dirac fermions at 
each KK level. This lead to the factor 2 appearing in C\ and C2 since the KK left 
and right- handed chiral states contribute to the closed fermion one loop diagrams. 
That is, 

GdUEDBulk = I Y20 51 ~8 92 T 53 J ^ ' ' ~ ' ' 

GuUEDBidk = (j^9l + y g\ + f flf) W) - 1) , (24) 

GeUEDBulk = + -§ 9 ^) ~ ^ ' 

auEDBuik = (S(t) - l)(4T t - -^gj - ^g\ + A 

GlUEDBulk = 2(>S(t) — 1) , 
G2JJ EDBulk — (S(t) — 1) . 
GeUEDBulk = — 1) ■ 

Following the convention of Ref jT5j. we use the coefficient (S(t) — 1) = — 1), 
which depends on the energy scale and also lead to a reproduction of the SM before 
crossing the threshold of the first KK state at [i — 1 /R. This coefficient corresponds 
to s = for the UED case only in Ref l43l .. 

Ho J 1 1 

5.3. UED SM Brane 

For the case where the fermions are restricted to the brane, we obtain the coefficients 
from Ref [33]. and Ref [43]. with a difference related to the SU(b) normalisation 
factor: 

GdUEDBrane = \ j9x + \sl + 2 (<S'(i) ~ 1) > 

GuUEDBrane = ( ^ff? + \& + 8ff|J 2(S(t) - 1) , (25) 

f 9 9 \ 

GeUEDBrane = [ |5l + 4 #2 J 2(S(i) - 1) > 

OyEUBranE = 2(S(t) - 1) f - 3#2 + A) , 
GlUEDBrane = > 
GiUEDBrane = 2(S'(t) — 1) , 
GdUEDBrane = 2(S(t) — 1) ■ 
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Note that the coefficient C\ — since we do not have a trace of fermionic loops as 
the fermions are restricted to the brane. 



6. Yukawa evolutions in the MSSM and 5D MSSM 

In the following we write the general form of the evolution equations for the various 
MSSMs we shall consider here, where we shall use a notation similar to the ones of 
Refs [18,45 .. Note that the beta functions contain terms quadratic in the cut-off, 
where this part dominates the evolution of the Yukawa couplings and of k. The 
top Yukawa coupling becomes non- perturbative before the gauge couplings thus 
limiting the range of validity of the effective theory. 

16tt 2 ^ = Y d \T d C -G d + (3YjY d + Y^Y u )c} , 

16^ 2 ^ = Y u \T U C - G u + (3y u tY„ + Y}Y d )c} , (26) 

16^ = ak + ([Y e T Y e *]k + k[YjY e }) C . 

where 

T d = 3Tr(Y^Y d )+Tr(Y^Y e ) , 

T u = 3 Tr(Y^Y u ) , (27) 
T e = 3 Tr(YjY d ) + Tr(Y}Y e ) . 



6.1. MSSM 

The MSSM, as a limiting case of the 5D models we shall consider in the following, 
and also when < t < ln( M x R ) (that is the energy we consider for the evolution 
from Mz to 1/R) the coefficients in the evolution equations are: 

rt I 2 i o 2 , 16 2 

GdMSSM — I -J^9l + 6 9l + -^93 

Gumssm = (^g\ + 3gi + ys 3 2 ) ! 2s 5 

GeMSSM — (^9l + 3.92 

62 2 
&MSSM = 2 T u - -g 1 - 6g 2 , 



5" 



Cmssm — 1 , 
Cmssm = 1 • 



These coefficients are modified when we enter the energy regime where the effects 
of the extra dimensions set in. The modifications depend on which particles are 
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decoupled and on the structure of the model. We shall consider two cases, one in 
which all particles can propagate in the extra dimensions (bulk case) and the other 
in which fermionic matter fields are constrained to the brane (brane case). 



6.2. Bulk 

When the energy scale E > 1/R or when the energy scale parameter t > hi( M * fl ), 
the coefficients in the 5D MSSM, for all three generations propagating in the bulk, 
can be expressed as: 



6.3. Brane 

However, when all matter superfields are constrained to live on the 4D brane, the 
quadratic evolution due to the sum over the two KK towers will be milder. We note 
also that the traces are the same in the MSSM because we don't have a fermionic 
loop in the extra- dimension in the brane case. The coefficients of the evolution 
equations are given by: 




(29) 



C^Dbulk 
C^Dbulk 



Ot5Dbulk 



2C 5Dbulk T u - (-gf + Qgl)S{t) , 
irS(t) 2 . 




(30) 




^bDbrane 



5Dbr (inc. 
'sDbrane 



2 T u - (^-gl + 9gl)S(t) , 
2S(t) , 



1 . 
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7. Scaling of the Yukawa couplings and the CKM matrix 

It is well known that in the SM, the quark sector's flavor mixing is parameterised 
by the CKM matrix: 

(V ud V us V ub \ /V n V 12 V 13 \ 
Vckm = V cd V cs V cb = V 21 V 22 V 23 , (31) 

\V td V ts VtJ \v 31 v 32 v 33 ) 

which makes it possible to explain all flavor changing weak decay processes and CP- 
violating phenomena to date, where the 10 year run of Babar at SLAC and the Belle 
detector at KEK has greatly improved our knowledge of the CKM matrix elements. 
In particular, for the standard parameterisation of the CKM matrix, which has the 
form: 

-S12C23 - Ci 2 S 2 3Sl3e I<5 C12C23 - Sl2S23Sl3e i<5 s 23 c 13 , (32) 

S12S23 - c 12 c 23 s 13 e lS -c 12 s 23 - s 12 c 23 s 13 e lS c 23 c 13 / 

where s\ 2 = sin #12, c\ 2 = cos #12 etc. are the sines and cosines of the three mixing 
angles Q\ 2l 9 23 and #13, and S is the CP violating phase. 

The CKM matrix arises from a consideration of the square of the quark Yukawa 
coupling matrices being diagonalised by using two unitary matrices U and V, 

diag(/ 2 ,/ 2 ,/ 2 ) =UY*Y«Ui , 

&mg(hlh 2 s ,hl)=VYlY d V^ , (33) 

in which / 2 , / 2 , / t 2 and h 2 d ,h 2 sl h 2 are the eigenvalues of Y^Y U and Y^Y d respectively. 
It follows that the CKM matrix appears as a result of the transition from the quark 
flavor eigenstates to the quark mass eigenstates upon this diagonalisation of the 
quark mass matrices: 

Vckm = UV^ . (34) 

From the full set of one-loop coupled RGE for the Yukawa couplings and the CKM 
matrix, together with those for the gauge coupling equations, one can obtain the 
renormalisation group flow of all observables related to up- and down-quark masses 
and the CKM matrix elements. 

The RGEs are very important tools to show the properties of the quark masses 
and the CKM matrix at different energy scales. We write down the general form 
for the evolution of f 2 , h 2 and the variation of each element of the CKM matrix 
V ik [H1E31I33] in the SM, the UED SM, the MSSM and the 5D MSSM. 
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7.1. SM, UED Bulk SM and UED Brane SM 

16.4 
dt 

„dh 2 
Wtt 2 —^- 
dt 



16^ 



dt 

,dVik 
dt 



fMTuA - G u ) + Wf 2 - 2Bj2h 2 \V tJ \ 2 } , 

3 

hp(T d A - G d ) + Wh 2 -2BJ2 fi\V*j\ 2 } , 



(35) 



16tH-^ = yi[2{T e A - G e ) + Wyj] , 



3 -B 
2 



h 2 4- h 2 



where A = 5 = 1 in the SM, A = 2S(t) - 1, B = S(t) in the UED Bulk SM and 
A = 0, B = 2S(t) in the UED Brane SM. The gauge couplings G for the SM, the 
UED Bulk SM and the UED Brane SM arc written in Eq.(|23]), Eq.([24]) and Eq. 
respectively. 



7.2. MSSM, 5D bulk and 5D brane 

>dff 



16?r 



dt 
dt 



16*^ 
dt 



= f 2 [2(T u C -G u ) + 6Cf 2 + 2CJ2 hjlVijf] 

3 

= h 2 [2(T d C - G d ) + 6Ch 2 + 2Cj2ff \Va\*] , 

i 

= y 2 [2(T e C-G e ) + 6Cy 2 } , 



16vr 



2 dVjk 
dt 



= C 



(36) 



f 2 4- f 2 h 2 4- h 2 

Eil_Jj_h 2 \r. v* t/ i l k m f 2 ir* 



P - f 



fl V rm V 3kV iT 



U2 J j jm v Jk Vlm 



where we use the same forms as in Ecis. (|28I29I30[) to fix the coefficients C, C and 
gauge couplings G to describe each model. 



8. Comparison of the models and implications 

In the following we compare the main results for the different models and their phys- 
ical meaning both in term of experimental limits and of the theoretical implications. 
However, as pointed out in the Appendix, there is a subtlety involved the running 
of the physical parameters in the UED bulk model. Although we plot their running 
up to the gauge unification scale, in fact, as illustrated in Figl27[ the introduction 
of new ultraviolet cutoff becomes imperative due to the scalar potential stability 
condition, and beyond this scale new physics should appear. For further discussion 
and quantitative analysis in the UED bulk model, refer to Ref [46] for details. In 
contrast, in the UED brane model, the physics parameters have a full running till 
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the gauge unification scale, since the Higgs self coupling evolution has a finite value 
which thus excludes the vacuum stability concern and validates the theory up its 
full scale. [33] 

8.1. Top Yukawa coupling 

UED SM: Bulk and Brane cases 

In Fig|4]the initial Yukawa couplings are given by the ratios of the fermion masses 
to the Higgs vacuum expectation value. The Yukawa couplings evolve in the usual 
logarithmic fashion when the energy is below 2 TeV, 8 TeV, and 15 TeV for the three 
different cases. However, once the first KK threshold is reached, the contributions 
from the KK states become more and more significant. The evolution of ft (see 
Eq. (|35|) ) depends explicitly on the cutoff A, which have finite one- loop corrections 
to the beta functions at each massive KK excitation level. Therefore, the running 
of the Yukawa couplings, or more precisely, the one-loop KK corrected effective 
four dimensional Yukawa couplings, begins to deviate from their normal orbits and 
start to evolve faster and faster. For the compactification radius = 2 TeV, the 
Yukawa couplings evolve faster than the other two, reaching its minimum value at 
the unification scale, after that point their evolution will "blow-up" due to the faster 
running of the gauge couplings and new physics would come into play. For the radius 
R^ 1 = 8 TeV, we find similar behaviour to the Br 1 = 2 TeV case, where the "blow- 
up" scale is not very far from that of 2 TeV case. However, for the third choice of 
radius, since the compactification radius is now much higher than the other two, we 
need more energy to push it further toward its "blow-up" point, which is at a higher 
unification scale. We also observe that the Yukawa couplings are quickly evolving 
to zero, however, a satisfactory unification of these seems to still be lacking. In the 
UED scenario, the unification of the Yukawa couplings is very desirable due to the 
fast power law running. As such, we have so far observed the Yukawa couplings 
all decrease with increasing energy, which agrees with what is observed in the SM, 
however, the Yukawa couplings are driven dramatically towards extremely weak 
values at a much faster rate. This is an interesting feature that distinguishes the 
UED model from that of the SM. 

5D MSSM Bulk 

The 4D MSSM contains the particle spectrum of a two-Higgs doublet model exten- 
sion of the SM and the corresponding supersymmetric partners. After the sponta- 
neous breaking of the electroweak symmetry, five physical Higgs particles are left in 
the spectrum. The two Higgs doublets H u and Hj, with opposite hypercharges, are 
responsible for the generation of the up-type and down-type quarks respectively. 
The vacuum expectation values of the neutral components of the two Higgs fields 

satisfy the relation v u 2 + vj 2 = f^f) = (174GeU) 2 . The fermion mass matrices 
appear after the spontaneous symmetry breaking from the fermion-Higgs Yukawa 
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couplings. As a result, the initial Yukawa couplings are given by the ratios of the 
fermion masses to the appropriate Higgs vacuum expectation values as follows: 

, rn UtC ,t , _ rn d , s ,b _ m e,/x,T 

Ju,c,t — ) "*d*8,b — ) V&.li.r — j I "J ' J 

V u V d V d 

where we define tan/3 = v u /vd, which is the ratio of vacuum expectation values of 
the two Higgs fields H u and Hd- 

From the complete sets of the RGEs we can run the renormalisation group flow of 
all observables related to up- and down-quark masses and the quark flavor mixings. 
For our numerical analysis we assume the fundamental scale is not far from the 
range of LHC scale, and set the compactification scale to be = 2 TeV, 8 TeV, 
and 15 TeV respectively. 

Actually, below the supersymmetric breaking scale the Yukawa and gauge cou- 
plings run in the usual logarithmic fashion, giving a rather slow change for their 
values. Therefore, for supersymmetric breaking theories around TeV scales, for sim- 
plicity, we take the supersymmetric breaking scale Msusy = Mz in the present 
numerical study, and run the RGEs from Mz up to the high energy scales for our 
three different compactification scales. 

Additionally, as illustrated in Fig. [SJ for the case of universal 5D MSSM, once 
the first KK threshold is crossed at fi = R^ 1 , the power law running of the various 
beta functions causes the Yukawa coupling to rapidly increase following the rapid 
increase in the gauge coupling constants in the left panel of FigJ3] From Eg . ([55)1 
and Eg. (|2T)1) we can find the quadratic term of S(t) providing a positive contribu- 
tion to the Yukawa beta functions, which is in contrast to beta functions of the 
gauge couplings (which include terms only linear in S(t)). Therefore, from Eg. p6[) . 
the positive contribution from S(t) terms will dominate the negative contributions 
from the gauge couplings, and cause the Yukawa couplings to increase rapidly. This 
behaviour can be observed for both small and large tan /3 cases. However, as il- 
lustrated in the first graph of Fig. [5j for small tan/3, the Yukawa coupling has a 




Fig. 4. The Yukawa coupling ft for the top quark in the UED SM as a function of the scale 
parameter t, for the bulk case (left panel) and the brane case (right panel) where the solid line is 
the SM evolution and for different compactification scales: _R — 1 = 2 TeV (red, dotted line), 8 TeV 
(blue, dot-dashed line), and 15 TeV (green, dashed line). 
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large initial value, therefore it blows up at a relatively low energy as compared with 
the case for large tan/3. As a result, as one evolves upward in the scale, the top 
Yukawa coupling is rising with a fast rate and is pushed up against the Landau pole 
where it becomes divergent and blows up. In the vicinity of this singular point the 
perturbative calculation becomes invalid, and the higher order corrections become 
significant. The Landau pole also indicates that there is an upper limit on the value 
of the gauge couplings where new physics must emerge before the Yukawa couplings 
diverge. 



5D MSSM Brane 

In the brane localised matter field scenario, the beta function has only linear terms 
in S(t), which is comparable with the S(t) term in the beta function for the gauge 
couplings. As depicted in Fig. [6j for a small value of tan/3, we have a large initial 
value of ft and the gauge coupling contribution to the Yukawa beta function is sub- 
dominant only. Therefore, as discussed previously, the Yukawa coupling f t increases 
rapidly as one crosses the KK threshold at n = resulting in a rapid approach 
of the singularity before the unification scale is reached. However, for an interme- 
diate value of tan /3, we have a relative smaller initial condition for the top Yukawa 
coupling and the Yukawa terms in the beta function become less important. The 
contributions from the gauge couplings may then become significant, which leads 
to a net negative contribution to the beta functions. Therefore, the curvature of 
the trajectory of the top Yukawa evolution might change direction, and the Yukawa 
evolution will decrease instead of increasing. This behaviour would become more 
obvious for a large value of tan j3. As observed in Fig. [6j for tan j3 = 30, we in- 
deed observe the decreasing behaviour of the top Yukawa couplings. This behaviour 
provides a very clear phenomenological signature, especially for scenarios with a 
larger tan (3 and that are valid up to the unification scale where the gauge couplings 



Runnings of f, for taii/5— 1 Runnings of f t for tan/?-30 



Fig. 5. The Yukawa coupling ft for the top quark in the bulk case of 5D MSSM as a function of 
the scale parameter t, for (left panel) tan /3 = 1 and (right panel) tan ft = 30 where the solid line 
is the MSSM evolution and for different compactification scales: ij -1 = 2 TeV (red, dotted line), 
8 TeV (blue, dot-dashed line), and 15 TeV (green, dashed line). 
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converge. 



8.2. CKM Matrix 

Because of the arbitrariness in choice of phases of the quark fields, the phases of 
individual matrix elements of the Vckm are not themselves directly observable. 
Among these we therefore use the absolute values of the matrix element | Vij | as 
the independent set of rephasing invariant variables. Of the nine elements of the 
CKM matrix, only four of them are independent, which is consistent with the four 
independent variables of the standard parameterisation of the CKM matrix. For 
definiteness we choose the \V u b\, \V c b\, \V US \ and the Jarlskog rephasing invariant 
parameter J — lmV u dV cs V* s V* d as the four independent parameters of Vckm- 



UED SM 

In Fig|7]we specially plot the evolution of \V u b\ for the UED bulk and brane cases. 
For the evolution of \V c b\ and |V^ S | we can observe similar behaviours, i.e., they 
all increase with the energy scale; the variation rate become faster once the KK 
threshold is passed. The absolute values of all the remaining magnitudes of the 
CKM matrix elements can be obtained from the unitarity equations, as depicted in 
FigUl with increasing energy the running of the CKM matrix shows a pronounced 
pattern at the point where the KK modes are excited, thus we could determine 
the renormalisation group evolutions of the full CKM matrix. As can be seen from 
Eq. ([35]) . the evolution of the CKM matrix is governed by the Yukawa couplings 
and the factor S{t). They evolve faster in the region where the power law scaling of 
the Yukawa couplings becomes substantial. Therefore, the renormalisation effect is 
explicit for mixings involving the third family, i.e., \V u b\ and \V c b\, due to the large 
value of their Yukawa couplings. Because of the smallness of the Yukawa coupling 



Runnings of f t for tan/3- 1 
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Runnings of f, for tan/J-30 




Fig. 6. The Yukawa coupling ft for the top quark in the brane case of 5D MSSM as a function of 
the scale parameter t, for (left panel) tan = 1 and (right panel) tan /3 = 30 where the solid line 
is the MSSM evolution and for different compactification scales: R~ 1 = 2 TeV (red, dotted line), 
8 TeV (blue, dot-dashed line), and 15 TeV (green, dashed line). 
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terms, the renormalisation group flow of the mixing between the first two families, 
i.e., the Cabibbo angle of \V US \, turns out to be very small. Although the mixing 
angle increases all the time, it is rather inert, even in the UED model. 

5D MSSM 

In Figsf5]anc![H]we plot the energy dependence of \V u b\ from the weak scale all the 
way up to the high energy scales for different values of compactification radii R~ l 
for the bulk 5D MSSM case brane localised matter fields case respectively. In these 
sets of pictures we consider two indicative choices of tan /?, that of tan ft = 1 and 
tan P = 30. 

The running of the CKM matrix is governed by the terms related to the Yukawa 
couplings. These Yukawa couplings are usually very small, except for the top Yukawa 
coupling (which could give a sizeable contribution). The CKM matrix element V u b — 
9\?,e~ %b can be used to observe the mixing angle, #13. It decreases with the energy 
scale in a similar manner regardless of whether tan j3 is small or large. However, for a 
large initial value of ft (small tan/3), the mixing angles have a more rapid evolution 
and end in the regime where the top Yukawa diverges and develops a singularity. 
Quantitatively we observe from these plots that the value of \V u b\ change by more 
than 50%. 

8.3. The Jarlskog parameter 

UED SM 

We next turn our attention to the quark flavor mixing matrix, especially the com- 
plex phase of the CKM matrix which characterises CP-violating phenomena. This 
phenomena has been unambiguously verified in a number of K — K and B — B 
systems. For the parameter J (Fig fTU|) . the characteristic parameter for the CP 
non-conservation effects, its variation becomes very significant. The larger the value 



().(K)4() 




t t 

Fig. 7. The CMK matrix elements \V u b\ in the UED SM as a function of the scale parameter t, 
for the bulk case (left panel) and the brane case (right panel) where the solid line is the SM, for 
different compactification scales: R~ 1 = 2 TeV (red, dotted line), 8 TeV (blue, dot-dashed line), 
and 15 TeV (green, dashed line). 
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of the compactification radius R, the faster J evolves to reach its maximum. We 
observe an approximate 30% increase for J at the unification scale compared with 
its initial value. 



5D MSSM 

From FigsfTTl in contrast, the Jarlskog parameter decreases quite rapidly once the 
initial KK threshold is passed. However, when tan/3 is large, we have a relatively 
longer distance between the initial and terminating energy track, the evolution of 
J can be driven towards zero or even further. Besides, as can be seen explicitly in 
Ref [32] ■ , the beta functions of the evolution equations of the CKM elements are up 
to the third order of the CKM elements, which are comparable smaller than that of 
Jarlskog parameter's quadratic dependence on the CKM elements. This fact then 



Runnings of |V u t,| for tan/5-1 




Runnings of |V u t,| for tan/?-30 




Fig. 8. The CMK matrix elements \V u b\ in the bulk case of 5D MSSM as a function of the scale 
parameter t, for (left panel) tan (3 = 1 and (right panel) tan /3 = 30 where the solid line is the 
MSSM evolution and for different compactification scales: = 2 TeV (red, dotted line), 8 TeV 

(blue, dot-dashed line), and 15 TeV (green, dashed line). 
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Fig. 9. The CMK matrix elements \V u b\ in the brane case of 5D MSSM as a function of the scale 
parameter t, for (left panel) tan (3 = 1 and (right panel) tan /3 = 30 where the solid line is the 
MSSM evolution and for different compactification scales: R~ 1 = 2 TeV (red, dotted line), 8 TeV 
(blue, dot-dashed line), and 15 TeV (green, dashed line). 
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leads to the relatively large variation of J with the increase of energy. Furthermore, 
for tan/3 = 30, the Jarlskog parameter drops almost to zero, which sets the effect of 
the SM CP violation to being very small. Note, however, that in a supersymmetric 
theory other sources of CP violation beyond the SM ones are typically present, 
therefore only a complete and detailed study of a specific model would allow us to 
establish the strength of the CP violating effects. 

For the matter fields constrained to the brane, in Figs l9l and fT2l we observe that 
the evolutions of these mixing angles and CP violation parameter are decreasing 
irrespective of whether the top Yukawa coupling grows or not. For small tan /3 we see 
similar evolution behaviours for these parameters as in the bulk case. The decreases 
in the value of V u b and J are much steeper, due to rapid growth of the top Yukawa 
coupling near the singular point. However, as tan/3 becomes larger, e.g. tan/3 = 30, 
the top Yukawa coupling evolves downward instead of upward. The decreases in 
these CKM parameters then becomes much milder towards the unification scale; 
though the reduction to effectively zero in the Jarlskog parameter persists. As a 
result, for the brane localised matter field scenario, it is more desirable to have a 
large tan/? for theories that are valid up to the gauge coupling unification scale. 

In summary, for the two 5D MSSM scenarios with matter fields in the bulk or on 
the brane, we have performed the numerical analysis of the evolution of the various 
parameters of the CKM matrix, and both cases give us a scenario with small or 
no quark flavor mixings at high energies, especially for the mixings with the heavy 
generation. The evolution equations which relate various observables at different 
energies, and also allow the study of their asymptotic behaviours, are particularly 
important in view of testing the evolution of the Yukawa couplings. In the universal 
5D MSSM model, the evolution of these CKM parameters have a rapid variation 
prior to reaching a cut-off scale where the top Yukawa coupling develops a singularity 
point and the model breaks down. For the brane localised matter fields model, we 
can only observe similar behaviours for small values of tan while for large tan /3, 
the initial top Yukawa coupling becomes smaller, the gauge couplings then play 





Fig. 10. The Jarlskog parameter J in the UED SM as a function of the scale parameter t, for the 
bulk case (left panel) and the brane case (right panel) where the solid line is the SM, for different 
compactification scales: R~ 1 = 2 TeV (red, dotted line), 8 TeV (blue, dot-dashed line), and 15 
TeV (green, dashed line). 
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a dominant role during the evolution of the Yukawa couplings, which cause the 
Yukawa couplings to decrease instead of increasing. As such the variations of these 
CKM parameters have a relatively milder behaviour, and the theory is valid up the 
gauge coupling unification scale. 



9. Neutrino parameter evolutions 

In a similar way to what done for quark parameters, we can study the evolution of 
the masses, mixing and phases in the neutrino sector. In the following we first estab- 
lish our conventions for these parameters and introduce the main present bounds 
and values from experiments. We than introduce the corresponding renormalisation 
evolution equations in the various models and discuss the numerical values obtained 
for the evolution of the parameters in different models, as previously done for the 
quark sector. 



Runnings of the CP violation J for tan/5- 1 




Runnings of the CP violation J for tan/5-30 




Fig. 11. The Jarlskog parameter J in the bulk case of 5D MSSM as a function of the scale 
parameter t, for (left panel) tan (5 = 1 and (right panel) tan f) = 30 where the solid line is the 
MSSM evolution and for different compactification scales: _R — 1 = 2 TeV (red, dotted line), 8 TeV 
(blue, dot-dashed line), and 15 TeV (green, dashed line). 
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Fig. 12. The Jarlskog parameter J in the brane case of 5D MSSM as a function of the scale 
parameter t, for (left panel) tan /3 = 1 and (right panel) tan f) = 30 where the solid line is the 
MSSM evolution and for different compactification scales: R~ 1 = 2 TeV (red, dotted line), 8 TeV 
(blue, dot-dashed line), and 15 TeV (green, dashed line). 
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9.1. Conventions for masses and mixing parameters 

The mixing matrix which relates gauge and mass eigenstates is defined to diagonalisc 
the neutrino mass matrix in the basis where the charged lepton mass matrix is 
diagonal. It is usually parameterised as follows [47] : 



C12C13 S12C13 si 3 e 

r = \ -S12C23 - c 12 s 2 3S 13 e~ tS c 12 c 23 - s 12 s 23 s 13 e lS s 23 c 13 
S12S23 - Ci 2 c 2 3Si3e z ' 5 -C12S23 - si 2 c 23 s 13 e lS c 23 ci 3 




with Cij = cosOij and Sy = smOij (ij = 12, 13,23). We follow the conventions of 
Ref. [45] to extract mixing parameters from the PMNS matrix. 

Experimental information on neutrino mixing parameters and masses is obtained 
mainly from oscillation experiments. In general Am 2 atm is assigned to a mass squared 
difference between v 3 and v 2 , whereas Am 2 sol to a mass squared difference between v 2 
and v\ . The current observational values are summarised in Table [2] Data indicates 
that Ai«j ol <C Am„ tm , but the masses themselves are not determined. In this work 
we have adopted the masses of the neutrinos at the Mz scale as m\ — 0.1 eV, 
m 2 = 0.100379 eV, and 7713 = 0.11183 eV, as the normal hierarchy (whilst any 
reference to an inverted hierarchy would refer to m 3 = 0.1 eV, with m 3 < mi < m 2 
and satisfying the above bounds) . For the purpose of illustration, we choose values 
for the angles and phases as the Mz scale as: 6*12 = 34°, #13 = 8.83°, $23 = 46°, 
5 = 30°, 0i = 80° and (j) 2 = 70°. 



Table 2. Present limits on neutrino 
masses and mixing parameters used in 
the text. Data is taken from Ref. I48| for 
sin 2 (20i 3 ), and from Ref. gj. 



Parameter Value (90% CL) 

sin 2 (2(?i 2 ) 0.861C™£) 

sin 2 (26*23 ) > 0.92 

sin 2 (20i 3 ) 0.092 ±0.017 

Am 2 aol (7.59 ± 0.21) x 10" 5 eV 2 

Am 2 t ° m (2.43 ± 0.13) x 10" 3 eV 2 



The evolution equation for the observables in our 5D MSSM are taken from [24] , 
As expected tan j3 plays an important role as all the mixing angles and phases 
depend on y T (see Appendix C in [23]). However, the new degrees of freedom (the 
extra-dimensional fields giving rise to KK excitations of the zero modes) become 
important at energies corresponding to their masses. In the following we study the 
evolution of the relevant parameters, such as Am^,, Am^ (m and the angles and 
phases, as a function of the energy scale and of tan (3. Only some selected plots will 
be shown and we will comment on the other similar cases not explicitly shown. 
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9 - 2 - A ™L and Am ltm 

For the UED SM, we see different behaviour for the brane case Fig[T3]and bulk case 
Fig. [T?] Once the KK threshold is reached, both Am^ and Am^ m decrease with 
increasing energy in the brane case, but they increase with the energy in the bulk 
case for the different radii of compactification. The evolution of masses depends on 
the evolution of y T and k coupling and the RG runnings in the UED SM bulk model 
are generally larger than those in UED SM brane model. This is due to the fact that 
the coefficient C\ = in the brane model (Eq.p5)l) and 2(S(t)~ 1) in the bulk model 
fEq. flMlO and also there is difference in a in the two equations due to the trace of 
charged-fermion Yukawa couplings in bulk model whereas such a contribution does 
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Fig. 13. Evolution of Am 2 sol (left panel) and Am^ im (right panel) as a function of the scale 
t = \n(fi/Mz) with matter fields constrained to the brane in the UED SM. The black line is the 
SM evolution, the red (small dashes) is for R' 1 ~ 1 TeV, the blue (dash-dotted) R~ l ~ 4 TeV, 
the green (large dashes) R' 1 ~ 15 TeV. 
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Fig. 14. Evolution of Am? o[ (left panel) and AmJ lm (right panel) as a function of the scale 
t = ln(fi/M z ) with matter fields in the bulk in the UED SM. The black line is the SM evolution, 
the red (small dashes) is for ~ 1 TeV, the blue (dash-dotted) R~ 1 ~ 4 TeV, and the green 

(large dashes) R _1 ~ 15 TeV. 
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not exist in brane model due to the absence of fermion KK excitations (see the T 
term in Eg . ((23|) ~) . This lead to the increasing of observables in the bulk case and the 
decreasing in the brane case. 

For the 5D MSSM, in general, in the brane case, the evolution has the same 
form for the three masses mi, rri2, W3. This leads to a reduction of up to a factor 
of two for the masses at t = 6 (for a large radius, R^ 1 = 1 TeV) with respect 
to the MSSM values at low energies (smaller radii give a weaker effect as the KK 
excitations contribute to the evolution equations at higher energies) . This prediction 
is extremely stable and can be explained as the evolution of the masses is governed 
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Fig. 15. Evolution of Am^ ol (left panel) and AmJ lm (right panel) as a function of the scale 
t = \n(p/Mz) with matter fields constrained to the brane for tan/3 = 30 in the 5D MSSM. The 
black line is the MSSM evolution, the red (small dashes) is for R~ 1 ~ 1 TeV, the blue (dash-dotted) 
R^ 1 ~ 4 TeV, the green (large dashes) R -1 ~ 15 TeV. 
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Fig. 16. Evolution of AmJ o| (left panel) and AmJ lm (right panel) as a function of the scale 
t = \n(n/M z ) with matter fields in the bulk for tan/3 = 30 in the 5D MSSM. The black line is the 
MSSM evolution, the red (small dashes) is for R" 1 ~ 1 TeV, the blue (dash-dotted) R' 1 ~ 4 TeV, 
and the green (large dashes) H _1 ~ 15 TeV. The evolution is towards a non-perturbative regime, 
where the Yukawa coupling develops a Landau pole and the effective theory becomes invalid. 
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by the equation 

1 



d m t (a + G iy 2 T ) , (38) 

where the coefficients Ci induce a non-universai behaviour and the parameter a 
contains the up Yukawa couplings and the gauge coupling terms (detailed in Ap- 
pendix C of Ref [24]. and Eq.([3"0)l), In contrast to the MSSM, the evolution in the 
brane case is completely dominated by the universal part. The essential point is 
that in the MSSM the positive contribution to a, approximately 6y 2 , is of the same 
order as the negative contribution from the gauge part. In our case the gauge part 
has a large pre-factor S(t) = e l MzR with respect to the MSSM which makes it 
completely dominant compared to any other contribution. As energy increases, we 
can write: 

1 



16tt 2 



<»'/r (|?i 2 + 9<7 2 2 ) M/i 



< . (39) 



From this approximation we immediately see that all masses decrease with increas- 
ing energy and eventually tend to zero if the evolution equations can be trusted up 
to a high energy. 

In the 5D MSSM bulk case the evolution has again the same form for the three 
masses mi, m^, but the behaviour is the opposite as the masses increase at 
high energy because all matter fields propagate in the bulk and contribute to the 
evolution. In detail this can be seen by the fact that even if the gauge part gets a 
large pre-factor S(t), the Yukawa part gets in this case a pre-factor S(t) 2 , which 
changes the sign of the derivative with respect to the previous case: 



16tt 2 



> . (40) 



From this approximation we see that all masses increase with increasing energy 
scale. 

The situation is more involved when analysing the mass squared differences. 
We plot in FigsfT5l and [TBI the evolution of Am 2 oi and Am 2 tm both for the matter 
fields on the brane and for all fields in the bulk for tan/3=30 and different radii of 
compactification. In the brane case different behaviours as a function of the energy 
scale are possible as a relatively large interval in energy range is allowed for the 
effective theory. As explicitly illustrated in Fig. [15l the relevant radiative corrections 
controlled by the gauge fields in Eq. (|39[) become dominant as energy goes up, which 
tends to reduce mass splitting, and an approximately degenerate neutrino masses 
spectrum at the high energy scale m\ ~ 7712 ~ "^3 becomes favourable. This is in 
contrast with the MSSM, where the neutrino mass splitting becomes large at an 
ultraviolet cut-off. Therefore, it is very appealing that the neutrino mass splitting at 
low energy could be attributed to radiative corrections resulting from a degenerate 
pattern at a high energy scale. In Figll6[ the bulk case tends to a non-perturbative 
regime, where the unitarity bounds of the effective theory are reached much faster 
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and only a much shorter running can be followed using the effective theory. As seen 
in Eq.([40|). the quadratic terms related to S(t) dominate during the fast evolution. 
As such, the neutrino mass splitting becomes even larger at a high energy scale. 

9.3. Mixing angles 

Concerning the evolution of the mixing angles, as can be seen in Figs fTTH^l in 
the UED bulk and brane cases, we have very small variation from the SM case 
because we don't have dependence on tan/3 and there is not a quadratic term of 
S(t) in the RGEs so there is no strong enhancement in the running. However the 
mixing angles variation is more significant in the 5D MSSM in which the largest 
effect is for 812, with changes of more than 70% possible for the brane localised 

Runnings of d t 1 for UED bulk model Runnings of H , ( for LIED brane model 
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Fig. 17. Evolution of #13 as a function of the parameter t = hx{fx/Mz) with matter fields in the 
bulk (left panel) and on the brane (right panel) respectively in the UED SM. The black line is the 
SM evolution, the red (small dashes) is for R' 1 ~ 1 TeV, the blue (dash-dotted) ~ 4 TeV, 

and the green (large dashes) R' 1 ~ 15 TeV. 




Fig. 18. Evolution of 613 as a function of the parameter t = Va(fj,/Mz) with matter fields in the 
bulk for tan j3 = 5 (left panel) and tan /3 = 30 (right panel) in the 5D MSSM respectively. The black 
line is the MSSM evolution, the red (small dashes) is for R^ 1 ~ 1 TeV, the blue (dash-dotted) 
R- 1 ~ 4 TeV, and the green (large dashes) R^ 1 ~ 15 TeV. 
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matter field scenario. As observed, due to the large quadratic term of S(t) in the 
beta function, the 9±2 has a rapid and steep variation in the bulk case. However, 
for the brane case, it has a relatively longer evolution track with the #12 then being 
pulled further down until the termination point (where the effective theory becomes 
invalid). In contrast, the running of #13 and 623 is much milder. As demonstrated 
in Figs. [18] and HH changes in the values of #13 vary only a couple of degrees. For a 
larger value of tan j3 we have a relatively large Yukawa coupling to t, which leads 
to a large magnitude for its beta function, resulting in a relatively large variation 
during the evolution. However, a running to #13 = cannot be observed in any 
situation. From the evolution behaviour of #13, one can see that the renormalisation 
group running effects or finite quantum corrections are almost impossible to generate 
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Fig. 19. Evolution of 613 as a function of the parameter t = ln(fi/Mz) with matter fields con- 
strained to the brane for tan = 5 ( left panel) and tan /3 = 30 (right panel) in the 5D MSSM 
respectively. The black line is the MSSM evolution, the red one (small dashes) is for ~ 1 TeV, 
the blue (dash-dotted) R^ 1 ~ 4 TeV, and the green (large dashes) ~ 15 TeV. 
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Fig. 20. Evolution of 612 as a function of the parameter t = hi(fj,/Mz) with matter fields in the 
bulk (left panel) and on the brane (right panel) respectively in the UED SM. The black line is the 
SM evolution, the red (small dashes) is for R^ 1 ~ 1 TeV, the blue (dash-dotted) i? -1 ~ 4 TeV, 
and the green (large dashes) R^ 1 ~ 15 TeV. 



33 



#13 = at a high energy scale, even though the power law enhanced evolution is 
considered during the running. Therefore, for the tri-bimaximal mixing pattern [49], 
in the current context with no other extreme conditions being taken into account, 
a slightly changed 6*13 could not be accommodated during the whole range of the 
energy scale. Similar trajectories are also observed for #23 (see Fig. [23]) . 

9.4. S phase 

The running of the Dirac phase 5 in the UED SM case is represented in Figl24l The 
variation is stable and similar for the bulk case (left panel) and brane case (right 
panel), there is very small deviation from the SM because all other mixing angles 
vary only by small quantities and the coefficient C which appear in the variation of 




Fig. 21. Evolution of 612 in the bulk (left panel) and on the brane (right panel) as a function of 
the scale t = ln(/x/Mz) for tan ,3 = 30 in the 5D MSSM. The black line is the MSSM evolution, 
the red one (small dashes) is for R~ 1 ~ 1 TeV, the blue (dash-dotted) R^ 1 ~ 4 TeV, and the 
(large dashes) ~ 15 TeV. 
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Fig. 22. Evolution of 623 as a function of the parameter t = hi(fj,/Mz) with matter fields in the 
bulk (left panel) and on the brane (right panel) respectively in the UED SM. The black line is the 
SM evolution, the red (small dashes) is for R^ 1 ~ 1 TeV, the blue (dash-dotted) i? -1 ~ 4 TeV, 
and the green (large dashes) R^ 1 ~ 15 TeV. 
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(5 (see Appendix D of [21]) are linear in S(t) and there is no dependence on tan/3. 

Noting that the Dirac phase 5 determines the strength of CP violation in 
neutrino oscillations. In the 5D MSSM, the runnings we include follow the gen- 
eral features presented in Figsj25] and [26l with large increases possible once the 
first KK threshold is crossed. From these studies we have seen that the vari- 
ation is bigger for high tan f3 with large changes appearing in the brane case 
when approaching the high energy scale. The recent results from the Daya bay 
and RENO reactor experiments have established a non zero values of $13. There- 
fore, the leptonic CP violation characterised by the Jarlskog invariant J ~ 
sin 9i2 cos O12 sin 623 cos 623 sin #i3COS 2 #i3 sin S becomes promising to be measured in 
the future long baseline neutrino oscillation experiments. As plotted, we can observe 
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Fig. 23. Evolution of 823 in the bulk (left panel) and on the brane (right panel) as a function 
of the parameter t = \n(/i/M z ) for tan fi = 30 in the 5D MSSM. The black line is the MSSM 
evolution, the red one (small dashes) is for R~ 1 ~ 1 TeV, the blue (dash-dotted) R~ 1 ~ 4 TeV, 
and the green (large dashes) R' 1 ~ 15 TeV. 



Runnings of S LiED bulk mode] Runnings of for ULD brauo model 

31.0 



10.5 
|. 10.0 
20.5 



:-" 50.0 



Fig. 24. Evolution of <5 as a function of the parameter t = ln(ii/Mz) with matter fields in the 
bulk (left panel) and on the brane (right panel) respectively in the UED SM. The black line is the 
SM evolution, the red (small dashes) is for R- 1 ~ 1 TeV, the blue (dash-dotted) R' 1 ~ 4 TeV, 
and the green (large dashes) R^ 1 ~ 15 TeV. 
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a relatively large evolution for the Dirac phase, even the maximum CP violation case 
<5 = ■§ could be achieved for relatively small input values. For leptogenesis related 
to the matter-antimatter asymmetry, we should note that the parameters entering 
the leptogenesis mechanism cannot be completely expressed in terms of low-energy 
neutrino mass parameters. Note that in some specific models the parameters of the 
PMNS matrix (which contains CP asymmetry effects) can be used [50I5T] . Here, the 
CP-violating effects induced by the renormalisation group corrections could lead to 
values of the CP asymmetries large enough for a successful leptogenesis, and the 
models predicting maximum leptonic CP violation, or where the CP-violating phase 
5 is not strongly suppressed, become especially appealing. Specific models with large 
extra dimensions in which leptogenesis is relevant at low scale can also be found in 
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Fig. 25. Evolution of the phase A as a function of the parameter t = \n(fi/Mz) with matter fields 
in the bulk for tan = 30 ( left panel) and tan /3 = 50 (right panel) in the 5D MSSM respectively. 
The black line is the MSSM evolution, the red one (small dashes) is for R~ 1 ~ 1 TeV, the blue 
(dash-dotted) R' 1 ~ 4 TeV, and the green (large dashes) R" 1 ~ 15 TeV. 
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Fig. 26. Evolution of the phase 8 as a function of the parameter t = ln(fi/Mz) with matter fields 
constrained to the brane for tan j3 = 30 ( left panel) and tan /3 = 50 (right panel) in the 5D MSSM 
respectively. The black line is the MSSM evolution, the red one (small dashes) is for R~ 1 ~ 1 TeV, 
the blue (dash-dotted) R' 1 ~ 4 TeV, and the green (large dashes) R 1 ~ 15 TeV. 
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Ref. [S2]. 

The running of the mixing angles are entangled with the CP-violating phases 
[24] . The phases <f>i and 02 do not affect directly the running of the masses, while the 
phase S has a direct effect on the size of dm/ dt, although its importance is somewhat 
reduced by the magnitude of #13. For further discussions of the correlation between 
these phases and mixing angles, refer to Refs [45l[53]. for details. 

Finally, whilst the above results and analysis were for the normal hierarchy of 
neutrino masses, we did also consider an inverted hierarchy, where from an anal- 
ysis of the equations presented in the Appendices of [24] we obtain the same fea- 
tures and results for neutrino mass runnings (though with different initial values 
at the Mz scale). As such, the figures for Am^ and Am^ tm remain unchanged. 
Possible changes in the angles and phases arise from the different signs for the 
(nij — mi)/(mj + m^ terms present in each evolution equation, where the 9 12 re- 
sults remain approximately the same (given the relative ordering of masses in these 
two hierarchies), and the small runnings of #13 and #23 would be up rather than 
down (though as already discussed, these runnings are quite small). 

10. Summary and Outlook 

The present review of the renormalisation group evolution of the masses, mixing 
angles and phases of the UED models in the quark and lepton sectors brings together 
the results we and other groups have obtained in the recent years in this subject 
using a common notation. The important physical points are discussed and the 
equations are written in compact way to show the unified approach to the different 
sectors of these models. For more technical details we refer to the existing literature. 

Concerning the UED standard model the evolution of the gauge couplings has 
a rapid variation in the presence of the KK modes and this leads to a much lower 
unification scale than in the SM. Due to the power law running of the Yukawa 
couplings, the rapid decrease of the Yukawa couplings with energy is in contrast 
to the logarithmic running predicted by the SM. The UED model has substantial 
effects on the hierarchy between the quark and lepton sectors and provides a very 
desirable scenario for grand unification. As for the energy scaling of the Jarlskog 
parameter J, the contribution of KK modes is substantial. Its numerical analysis 
shows that its variation can be raised to more than 30%. The scale deviation of 
renormalisation curves from the usual SM one depends closely on the value of the 
compactihed radius R. The smaller the radius is, the higher the energy scale we need 
to differentiate the UED curve from the SM one. A comparison between theoretical 
predictions and experimental measurements will be available once the LHC will be 
running at its full centre of mass energy. This will set limits on the parameters of 
the UED model, and a precise determination of J, \V u b\ or \V c b\ at high energy may 
lead to a discrimination between the SM and extra dimensional models. 

In the case of the 5D MSSM, we have reviewed the behaviour of the evolution 
equations for the quark and neutrino sector in a minimal supersymmetric model 
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with one extra-dimension. For quarks, the 5D MSSM scenarios with matter fields in 
the bulk or on the brane, give both results with small or no quark flavor mixings at 
high energies, especially for the mixings with the heavy generation. The evolution 
equations are particularly important in view of testing the evolution of the Yukawa 
couplings. The evolution of these CKM parameters have a rapid variation prior 
to reaching a cut-off scale where the top Yukawa coupling develops a singularity 
point and the model breaks down. For the brane localised matter fields model, we 
can only observe similar behaviour for small values of tan j3, while for large tan /3, 
the initial top Yukawa coupling becomes smaller, the gauge couplings then play 
a dominant role during the evolution of the Yukawa couplings, and therefore the 
Yukawa couplings decrease instead of increasing. As such the variations of these 
CKM parameters have a relatively milder behaviour, and the theory is valid up 
the gauge coupling unification scale. Concerning the neutrino sector, the evolution 
equations for the mixing angles, phases, Am^, and Am^ tm , within the two distinct 
scenarios, is also considered. A larger tan (3 typically leads to larger renormalisation 
group corrections. Neutrino masses evolve differently in the two models due to the 
sign of the (different) dominant contributions in the bulk and in the brane cases. 
For the brane case we find the approximate degenerate neutrino mass spectrum 
becomes more favourable at the ultraviolet cut-off. In the bulk case, the neutrino 
splitting becomes even more severe as the unitarity bounds of the effective theory 
are reached faster. Contrary to the large renormalisation effect of 6*12, the runnings 
of and 023 are relatively mild. We found a non-zero value for $13 during the 
evolution, which has no appreciable RGE running effects, even when power law 
evolution effects are considered. Therefore it is necessary to introduce new physics 
effects in order to achieve the tri-bimaximal pattern. The maximum CP violation 
case, 5 = |, could be achieved starting from a relatively small initial value. In 
general we can see that radiative effects may have a significant impact on neutrino 
physics. A non-zero Jarlskog invariant, which measures the magnitude of leptonic 
CP violation( expected to be measured in future long baseline neutrino oscillation 
experiments), could open the door for measurable CP violation in the leptonic 
sector. 
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Appendix A. 

In the SM, the Higgs boson mass is given by mjj = \f\v, where A is the Higgs 
self-coupling parameter and v is the vacuum expectation value of the Higgs field 
(where v = (V^Gp)^ 1 ^ 2 = 246 GeV is fixed by the Fermi coupling Gp). From 
the requirement that the scalar potential energy of the vacuum be bounded from 
below, the quartic coupling A should be positive at any energy scale. If ran is too 
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small, A becomes negative at certain energy scales and then induces a false and 
deep minimum at large field values, destabilising the EW vacuum. Therefore, above 
that scale, the validity of the SM is expected to fail and must be embedded in some 
more general theories that give rise to a wealth of new physics phenomena. 

In the five dimensional UED model compactified on a circle of radius R with a 
Z2 orbifolding, motivated by the new bounds on the mass of the SM Higgs boson 
around 125 GeV of ATLAS and CMS collaborations [54], we quantitatively analyse 
the Higgs self coupling evolution from the EW scale up to the unification scale and 
exploit its evolution behaviours for different compactification radii. 

The kinetic term of the scalar doublet has the following forms in the 5D UED 
model: 

TlR 

C Higgs = J dy(D M <P(x,y)yD M <P(x,y) , (A.l) 


where the covariant derivative is 

D M $(x, y) = I d M + ig\T a W M + ^?%) ®(x, y) , (A.2) 

where the gauge fields G M (x,y) f G M^-^> Wu{x,y) and B M (x,y) 

refer to the SU(3), SU(2) and U(l) gauge groups respectively. Note also the five 
dimensional gauge coupling constants <?|, g\ and g\ are related to the four dimen- 

g 5 

sional SM coupling constants (up to a normalisation factor) by gt = 1 , and 

V7T-R 

A 5 

similarly for the quartic coupling A 



I-kR 

The evolution of the Higgs quartic coupling is given by the beta function as 
follows 

= + ^UED ^ (A _ 3) 



where 



dt 



- 4Tr[3{Y^Yu) 2 + 3{Y^Y D ) 2 + (Y^Y E ) 2 } (A.4) 
where T t = Tr[3Y^Y v + 3Y^Y D + Y E Y E ] and 

= (S(t) 1) |l2A 2 - 3 + 3 52 2 ) A + QLrf + \g\g\ + igi 

+ 2(S(t) - 1) {4A T t - 4Tr[3(yjF a ) 2 + l^Yof + (Y^Y E ) 2 }} (A.5) 

for the UED bulk model, [T0I[5B] where both the fermion fields and the boson fields 
can propagate in the bulk, and 
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Fig. 27. The evolution of the Higgs self coupling for mu = 125 GeV, where the solid (black) line 
is the SM, the upward [downward] dotted (red) line is the R' 1 = 2 TeV UED brane model [UED 
bulk model], the upward [downward] dotted-dashcd (blue) line is the H _1 = 8 TeV UED brane 
model [UED bulk model], and the upward [downward] dashed (green) line is the R~ 1 = 15 TeV 
UED brane model [UED bulk model] . The ultraviolet cutoff incurred from the vacuum instability 
are observed around t = 4.5, 5.7, 6.3 for = 2, 8, 15 TeV respectively in the UED bulk model. 

for the UED brane model, [33] where the SM matter fields are localised to the brane, 
whilst the gauge and Higgs fields are propagating in the bulk. Note that the pure 
gauge terms in Eq. (|A.6j) arc different from these in Ref. [43, 57 . Explicitly, these 
terms cannot exactly resemble what is in the SM, since there are extra contributions 
from the couplings of the Higgs field with the A 5 component in the bulk space as 
well. 

In Figl27l starting from the initial value of run = 125 GeV, we focus on the 
evolution of the Higgs self-coupling and explore its behaviour and constraints on the 
compactification radii for which the validity of the theory is satisfied. For interesting 
values of a Higgs mass around 125 GeV, we can observe that, in the whole range 
from the EW scale up to the gauge unification scale, in the UED brane model 
the Higgs self coupling \(t) remains positive and its trajectory goes upward but 
remains finite when approaching the unification scale, whilst for the SM and the 
UED bulk model X(t) evolves towards a zero value before reaching the unification 
scale, which then incurs the vacuum instability and introduces an ultraviolet cutoff 
for the theory. 
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